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We present a model for fully developed triplet proximity effect in superconductor-ferromagnet 
heterostructures. Within the circuit-theory approximation, we evaluate the Green's functions, the 
density of states, and the Josephson current that depend essentially on the magnetization configu- 
ration. 
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Research of heterostructures that combine supercon- 
ducting (S) and ferromagnetic (F) elements gives insight 
to the problem of mutual influence of superconductiv- 
ity and ferromagnetism, allows realization of exotic su- 
perconducting states such as the L-O-F-F state and 
triplet ordering, and promises applications that uti- 
lize the spin degree of freedom. While this research 
has started more than three decades ago, interest to 
the above topics has yet stipulated recent important de- 
velopments, both theoretical and experimental. Those 
concern Josephson 7r-junctions, triplet supercon- 
ductivity, 0,0,0 an d Josephson spin valves. 00 

The superconducting proximity effect is characteris- 
tic for most S/F heterostructures and is distinct from 
the effect in non-magnetic S/N systems. In the normal 
part of a S/N structure, the superconducting correlations 
persist at distances of the order of normal-metal coher- 
ence length £jy This length scale can grow large at suf- 
ficiently small temperatures T. In a diffusive material, 
£at = \/hD/2nkBT, D being the diffusivity. In con- 
trast to this, superconducting correlations in a ferromag- 
net, where an exchange field h is present, are quenched 
at much shorter scale £/, = y/hD/h. Hence one might 
conclude that the proximity and Josephson effects are 
strongly suppressed in S/F heterostructures. However, 
some experiments |ldj seem to contradict this statement, 
indicating proximity correlations at much larger scale. 
Though these experiments may be explained by inter- 
face effects, they have motivated a proposal of an 
interesting mechanism for long-range proximity effect in 
ferromagnets 0-0] It was shown that inhomogeneity in 
the direction of exchange field generates superconduct- 
ing correlations of two electrons with the same spin, i.e. 
triplet correlations. Such triplet proximity effect (TPE) 
is not suppressed by an exchange field and penetrates the 
ferromagnet at the scale of fjy Recently, a substantial 
Josephson current has been reported for a fully polarized 
ferromagnet. |j| The experiment can only be explained 
by TPE. 

An immediate problem is that the theoretical predic- 
tions so far have been elaborated in assumption that TPE 
is weak and can be treated perturbatively. This makes it 
difficult to determine an unambiguous experimental sig- 
nature of TPE to distinguish it from the conventional ef- 



fect. Experimentally, the Josephson current due to TPE 
7] does not seem to be smaller than that due to a fully 
developed conventional proximity effect. 

In this work, we address a fully developed TPE that 
significantly changes the density of states (DOS) at the 
Fermi level. We show that the DOS increases. This is 
in contrast with complete suppression of DOS by a fully 
developed conventional proximity effect. Similar to the 
conventional effect, the change in the DOS is restricted 
to the energy window ~ Exh, the Thouless energy of 
the structure, provided Etk <C A, A being the energy 
gap in the superconductor. Therefore the corresponding 
Josephson current is of the same order as for the conven- 
tional proximity effect, though its magnitude essentially 
depends on the magnetization configuration in the struc- 
ture. For magnets where both spin directions are present 
at the Fermi surface (ferromagnetic metals) , we find both 
7r-junctions 00 (with negative supercurrents) and com- 
mon 0-junctions (with positive supercurrents), depending 
on the magnetization directions. For fully polarized mag- 
netic materials (half- metals) , we find a continuous de- 
pendence of the equilibrium superconducting phase dif- 
ference on magnetization directions. Finite current may 
be induced thereby at fixed zero phase difference. 

We concentrate on a S/F/S heterostructure fabricated 
by deposition of two superconducting electrodes onto 
a conducting ferromagnetic film (Fig. l a) .Ml We pro- 
ceed with the so-called circuit theory jl^l that is a 
finite-element technique for semiclassical Green's func- 
tion method L3J, which has been applied to S/F struc- 
tures in Ref. fcj. Circuit theory is convenient since we 
aim at presenting an idealized TPE without considera- 
tion and subsequent optimization of concrete geometry 
of the structure. Besides, it allows simple analytical pre- 
sentation of the results. 

We build a minimal circuit-theory model for a ferro- 
magnetic metal, later adjusting it to a half-metal. The 
left and right parts of the structure contain regions where 
superconducting and magnetic correlations meet (points 
1 and 3 in Fig l.a). Following 0, we represent each 
region by a normal-metal node connected to a supercon- 
ducting reservoir and a ferromagnetic insulating reser- 
voir (FIR). The role of FIR is to represent the exchange 
field hk (k = 1, 2, 3) induced in the node. The middle of 



2 



1 



a) y J 



F * : 



V -3 



b) 



4 





1 




1 



IfF 



FIX 



c) 



m A 



n, 



S/F) 



FIG. 1: Building up the circuit-theory model, a) S/F/S struc- 
ture exhibiting TPE Q- b) The corresponding circuit consists 
of three nodes representing the parts 1 — 3 of the structure, 
c) The nodes 1 and 3 are replaced by effective S/F reservoirs. 



the structure (point 2) is represented as node 2 connected 
to the nodes 1,3, and another FIR. To enable TPE, we 
allow arbitrary magnetization directions of all FIR's. It 
is assumed for simplicity that all connectors are of tunnel 
nature. The finite volume of each node and related deco- 
herence between electrons and holes is taken into account 
by introducing "leakage" matrix currents |l^ inversely 
proportional to the mean level spacing Si in each node. 
This defines a circuit presented in Fig. lb. 

The relevant variables are (retarded) Green's func- 
tions 01,2,3 at the nodes. These Green's functions are 
energy-dependent 4x4 matrices in the Nambu and spin 
space. They are determined from the equations that ex- 
press conservation of the "matrix currents" at each node. 
These equations read (k — 1,3) 
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Here, Gq = e 2 /irh is the conductance quantum, r and 
<x are vectors composed of Pauli matrices in the Nambu 
and spin space, respectively. The first term in these 
equations is the matrix current into the corresponding 
FIRs. The magnetization direction in FIR k is given by 
a unit vector n^. G s k is a Green's function of supercon- 
ducting reservoir k. We assume for simplicity that the 
relevant energies are much smaller than the supercon- 



ducting gap in the reservoirs, so we can use the energy- 
independent G s k = T\ cos ifik + T 2 sin ipk , ifik being the su- 
perconducting phase. Each superconducting reservoir k 
is connected to node A; by a tunnel junction with con- 
ductance gl . It would induce a proximity effect mini-gap 
Afc = gtSk/2GQ in node k if no other connections were 
present. The parameters g" 3 give the conductances of 
connections between the normal nodes. 

The two energy scales of the model are the typical 
exchange field h and the Thouless energy Eth — 9 n 5/GQ. 
The condition h < Etk corresponds to a short S/F/S 
structure with dimensions shorter or of the order of £/j . In 
this case, the conventional proximity effect overshadows 
TPE: singlet anomalous components of G in all nodes 
either exceed or are of the same order as the triplet ones. 
Since we want to single out the TPE, we turn to the 
opposite limit Etk *C h of longer structures. 

The best separation between the islands is achieved in 
the limit S 2 3> £1,3, where the regions 1 and 3 adjacent 
to the superconductors are much longer than the middle 
part of the system. We will see that in this limit the 
superconducting correlations in the nodes 1 and 3 are 
between electrons of opposite spin. As to the node 2, the 
correlations are between electrons of the same spin only: 
there TPE is present in its purest form. 

Under these conditions, the last term in the first 
Eq. can be neglected and G\^ can be determined sep- 
arately from G 2 ■ Owing to their big size, the nodes 1 and 
3 act as effective reservoirs for the node 2 (Fig. lc). The 
Green's functions at these nodes are determined by the 
competition between the corresponding superconducting 
reservoir and FIR and read 0] (fc = 1, 3; hu > A&; t± — 
[n±iT 2 }/2): 



G k 



h k r 3 - i(A* k T + + A fc T_)(ra fc • cr) 

VH-\A k \ 2 



(2) 



where we assume e <C hk,A.k- As seen from the struc- 
ture of |J2J|, the effective S/F reservoirs supply supercon- 
ducting correlations that are different for opposite spin 
directions. These correlations are most pronounced if 
hk — Afc. We will see that this condition is optimal for 
TPE. At this stage of research it is difficult to immedi- 
ately relate this condition to specifics of the structures in 
hand such as geometry, film thickness, transparency of 
S/F interface etc. Since the observations of Q suggest 
that TPE is close to the optimal one, we are convinced 
that this condition is realizable. 

Now the node 2 is connected to reservoirs only. Its 
Green's function is determined by the balance of the ma- 
trix currents into these reservoirs. The two S/F reser- 
voirs are connected in parallel, so that their net effect 
is additive and can be represented by a matrix M — 
(<7?Gi + 53"G 3 )/(5l 1 + S3 1 ) = M t 3 - i(M*r + + Mr_) • a. 
It is important that the resulting G 2 splits into two in- 
dependent blocks corresponding to two spin projections 
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on n.2- The separation into blocks allows us to treat 
half-metals on equal footing with ferromagnetic metals. 
While for a ferromagnetic metal both blocks contribute 
to physical quantities, only a single block does so in a 
half- metal. It is implied that in any case the S/F reser- 
voirs support both spin directions, 0] otherwise they 
would not be superconducting. 

We write the block structure as follows: 
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S = (1 - r 3 )/2 - ia y (n 2 ■ a)(l + r 3 )/2. 



In these notations, 
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-iM y \ 2 and the z axis 



where Af = (M -ie/E Th ) 2 - \M X 
is chosen in the direction of «2 ■ The result for G^ , j4j is 
obtained by replacing M^, ± iM 9 with M x =F iM y . The 
advantage of such notations is that the block structure 
is made explicit. The non-diagonal elements of G| cor- 
respond to triplet anomalous averages (t/'tV'tK (V'{V'{)> 
whereas the singlet-pairing averages V'i ) j van ~ 
ish. This manifests a pure TPE. The fully developed 
effect is characterized by hk |Afc|, in which case the 
non-diagonal and diagonal elements of Gfj are of the 
same order of magnitude. 

The triplet anomalous averages (ipiipi) ((V'j.V'l)) ac ~ 
quire a phase factor e* x (e~~ lx ) upon rotation by an an- 
gle x about the z axis. This leads to an interesting 
interplay between the superconducting phase difference 
03 — (f>i = <fi and the relative longitude angle between ri\ 
and n 3 , X3 ~ Xi = X- To see this explicitly, we express 
the magnetization direction vectors rii 3 in the spherical 
coordinates (x, standing for the longitude and latitude 
respectively). In these notations, 
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where a k = gi \ A k \ sin 9 k / (g? + g% ) yfh\ - |A fc | 2 , so that 
the superconducting phase and the longitude always 
come together. 

The DOS is determined from Gi and generally is dif- 
ferent for the opposite spin directions, 
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where ^q'^ is the normal-state DOS. As shown in 
Fig. Ufa), v always exceeds the normal-state value at 
£ "C Eth- This is in contrast with a suppression of v 
manifesting the conventional proximity effect. Such an 
enhanced DOS is therefore a signature of TPE. The peak 
at small energies is followed by a wider dip at energies 
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FIG. 2: The DOS in the central node a) versus the energy 
for = and \ = where 6\ = 63 increases from for the 
lower curve to tt/2 for the upper one with a step tt/8. b) The 
total DOS for a ferromagnetic metal at zero energy versus <j> 
for 9i — 6z — tt/2 where \ increases from to ■k with a step 
7r/4. Upon increasing \ a single peak at cf> = splits in two. 
The thick curve corresponds to a half-metal with x — k/2. 
For both graphs, /ii, 3 /Ai, 3 = 18/17. 



of the order of the Thouless energy Etk so that the to- 
tal number of states remains unchanged [Fig. Ufa)]. The 
dependence of v on the superconducting phase/longitude 
is also characteristic of TPE. As seen from Eq. JSJl, for a 
given spin direction the effect of superconducting phase 
can be always compensated by a rotation of one of the 
magnetization directions ni i3 about While for a fer- 
romagnetic metal the total DOS v = v> + is an even 
function of (j> and x, [Fig. Hfb)] this is generally not so 
for a half-metal. In the latter case, for a symmetric setup 
ai = a 3 , the DOS i/(0) can be modulated by either phase 
from 1 to its maximum value 1 / \/l — 4af /M5 2 . 

In the model under consideration, the Josephson cur- 
rent is given by the Keldysh component of the matrix 
current through either junction g™ 3 integrated over en- 
ergy, 
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FIG. 3: The Josephson current for a ferromagnetic metal ver- 
sus (j) for different x- a ) varying from to 7r with a step n/12. 
The first harmonic dominates the current except the vicinity 
of x = k/2 where it changes the sign, b) The vicinity of "0"- 
'V transition. Here x varies from tv/2 — n/36 ('V'-state, 
dotted line) through tt/2 - tt/72, tt/2, tt/2 + vr/72 (bistable, 
thin lines) to n/2 + n/36 ("0" -state, thick line). 
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The block structure of Gi gives rise to two contributions 
to the current that correspond to opposite spin direc- 
tions. The integral in Eq. J7J is logarithmic, converging 
at energies E Th < e < A cut , A cut ~ min(A, h,h- A). 
We assume T <C £r/i- With logarithmic accuracy, 
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Both the scale of the current (/ ~ gExh/e) and the log- 
arithmic structure are similar to the common proximity 
effect. The difference is the dependence of the current on 
the magnetization directions. For a ferromagnetic metal, 

Interestingly, the sign 



the total current / 



-I± is odd in < 



of the Josephson current is opposite to that in a common 
Josephson contact provided xl < it/ 2. This signals a 
7r-j unction [jUilllSj. which can be changed to a common 
0-junction by changing \. Since the accuracy of the log- 
arithmic approximation is always questionable, we eval- 
uate the integral numerically assuming A cut /ETh = 50. 
The resulting dependence of the current on both phases 
is not precisely harmonic, though close to it. (Figs. HI 111 . 
As seen in Fig. E^b), the second harmonic of the cur- 
rent becomes dominant in the vicinity of x = 7r/2. This 
implies that the transition between "ir" and "0" states 
follows a scenario of Ref. where both states are stable 
in this vicinity. 

For a half- metal, the situation is very different. In 
this case, the superconducting phase corresponding to 
the energy minimum just follows x arL d changes continu- 
ously instead of "jumping" between the values of and 
7r. There is a finite supercurrent at zero phase difference. 
Therefore, rotation of the magnetizations 711,3 about 
is equivalent to the effect of extra magnetic flux x^o/27r 
enclosed in a large loop that includes the junction. In 
our opinion, this facilitates an unambiguous experimen- 
tal verification of TPE. 

In both cases, not only the supercurrent is zero at the 
energetically favorable phase difference, but also TPE is 
reduced, vanishing completely for an ideally symmetric 
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FIG. 4: Tne Josephson current for a half-metal versus <f> for 
X varying from to tt/2 with a step n/8. The I — <f> curves 
shift to the left with increasing x- 



setup. This is because TPE increases the electron energy 
in the central node, as seen from the fact that the DOS 
is enhanced at low energies. 

In conclusion, we have proposed a simple model for 
fully developed triplet proximity effect in S/F/S struc- 
tures. In contrast to the common proximity effect, TPE 
enhances the DOS at low energies. The Josephson cur- 
rent exhibits a peculiar dependence on the magnetization 
configuration that is essentially different for a ferromag- 
netic metal and a half-metal. Those are signatures of 
TPE to be observed experimentally. 
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